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Abstract 

Credit Suisse First Boston (CSFB) launched in 1997 the model CreditRisk + which aims 
at calculating the loss distribution of a credit portfolio on the basis of a methodology from 
actuarial mathematics. Knowing the loss distribution, it is possible to determine quantile- 
based values-at-risk (VaRs) for the portfolio. An open question is how to attribute fair VaR 
contributions to the credits or loans forming the portfolio. One approach is to define the con- 
tributions as certain conditional expectations. We develop an algorithm for the calculations 
involved in this approach. This algorithm can be adapted for computing the contributions 
to the portfolio Expected Shortfall (ES). 



1 Introduction 



The CreditRisk + model ( |CSFB, 1997 ) has found wide-spread applications for the measurement 



of risk in credit portfolios. Reasons for the success of the model might be its free availability, the 
speed of calculations based on it due to fact that no Monte Carlo simulations need to be per- 
formed, and its detailed documentation. The basic idea in the model is to apply a methodology 
from actuarial mathematics. 

As soon as the loss distribution has been computed, it is an easy task to calculate the portfolio 
value-at-risk (VaR) . A further step might be to perform a risk diagnostics in the spirit of Litter-| 



man (199fj| ) . Such a diagnostics presupposes that the total portfolio risk can be attributed to the 



portfolio components in a way that detects the impact on the portfolio risk by the components. 
A first suggestion for determining such VaR contributions was made in the CreditRisk + doc- 



umentation (based on a decomposition of the portfolio variance (see |CSFB, 1997|) ; cf. also 



Biirgisscr et al., 1999) ). Further proposals can be found in Lehrbass et al. (2001 , sec. 5) or 
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Martin et al. (2001 ). In both these latter proposals the VaR contributions are motivated as ap- 
proximations to partial derivatives of VaR. Here we follow an alternative approach which relies 
on the fact that under certain continuity assumptions the partial derivatives of VaR with respect 
to the weights of the assets are just conditional expectations. This leads to an alternative defi- 
nition of VaR contributions (see (|2.7| )) which might be useful for checking the results according 
to the approaches in Lehrbass et al. (2001 ) or Martin et al. (200l| ). Moreover, this definition can 
easily be adapted for the case of Expected Shortfall (ES) (see ( f4.2| )) as risk measure. 



This text is based on the representations of the CreditRisk + model as given by Lehrbass et al 
(200l|) and |Gordy (200T|) . In section |, we introduce some basic features of the CreditRisk + 



model and provide motivation for the kind of VaR contributions to be discussed below. Then, 
in section El, we study the model in a more detailed manner in order to derive the main result 



of this paper, Corollary |3-4 It states that once it is possible to calculate the probability masses 
of the loss distributions, the computation of VaR contributions can be performed essentially 
the same way. This result is adapted to the case of ES in section f|. We conclude with a short 
summary of the results. 



2 Abstract description of the problem 

The key idea in CreditRisk + . Consider a portfolio with n loans or credits. Default or non- 
default of loan i, i = 1, . . . , n, by a fixed time horizon is indicated by the random variable D{ 
which can only take the values 1 for default or for non-default. Denote by V{ the exposure of 
loan i, measured as a positive integer multiple of a certain currency unit, pi = P[-Dj = 1] € (0, 1) 
is the default probability of loan i. The realized total loss Lq of the portfolio is then given as 

n 

L = ^ViDi. (2.1) 

i=i 

The following assumption is quite common for the purpose of dependency modeling in a portfolio 



of possible dependent assets (see Frey and McNeil, 2001 ). Suppose that all the random variables 



under consideration are defined on the underlying probability space (f2,.F, P). 

Assumption 2.1 There is a a-algebra A C T such that the default indicators D±, . . . ,D n are 

independent conditional on A. 

We will see examples for the choice of A in section ||. A may be regarded as the collection of 
systematic factors influencing the portfolio. 

If the default probabilities pi are small, the total loss distribution might not change too much 



when the Bernoulli variables D{ in ( 2A ) are replaced by random variables N\ , . . . , N n with 
non-negative integer values that conditional on A are independent and Poisson-distributed with 
A- measurable intensities R\ , . . . , R n > such that 

E[Ri] = Pi , i = 1,... ,n. (2.2) 
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In the CreditRisk + model, one does not compute the distribution of the random variable Lq 
from ( |2.lD but the distribution of the random variable 

n 

L = ^ViNt (2.3) 

4=1 

which, of course, is only an approximation to the true distribution. 



Which risk contributions should be considered? For a fixed (close to 1) level 5 E (0, 1), 

p{L) = q S (L) = inf{/ el: P[L < 1} > 6} (2.4) 

is a widespread measure of portfolio loss risk. It is called value-at-risk (VaR) at level 5 of L. In 
order to identify sources of particularly high risk in the portfolio, it is interesting to decompose 
p(L) into a sum Yli=i Pi^fi) wnere the risk contributions Pi(L) should correspond in some sense 
to the single losses Vi iVj. 



Litterman "(1996) suggested the decomposition 



i=l 



(2.5) 



h=0 



Decomposition (2.5) holds whenever the risk measure p is positively homogeneous (i.e. p(hL) = 
h p(L) for h > 0) and differentiable. Unfortunately, the distribution of the portfolio loss L, 
specified by (|2.3j), is purely discontinuous. Therefore, the derivatives of q$(L) in the sense of 



(2.5) will either not exist or be 0. 

Nonetheless, for a large portfolio even the positive probabilities P[L = I] will be rather small. 
On an appropriate scale, the distribution of L will therefore appear "almost" continuous. Taschej 



(19991) showed 

• that [Littermaril 's proposal Q2.5| ) is the only decomposition which is compatible with RO- 
RAC (Return on risk-adjusted capital) based portfolio management, and 

• that in case p(L) = qs(L) under certain continuity assumptions on the distribution of the 
asset losses one would have 



dh 



(h Vi Ni + L) 



v l E[N i \L = q 5 (L)}, 



1, ... ,n. 



(2.6) 



h=0 



Note that under the assumptions of Tasche (1999[) the event {L = qs(L)} has probability 0. 
Hence, the conditional expectation in (|2.6| ) has to be understood in the non-elementary sense 
(see purrett, 1996| , ch. 4). 
I I . . 

Litterman's and |Taschej 's considerations suggest the definition of 

Pi (L) = u i E[N l \L = q 5 (L)], i = l,...,n, (2.7) 
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as VaR contributions of the loans i = 1, . . . , n in the portfolio. Since P[L = 95 (L)] is positive by 
definition (^J) of VaR, the conditional expectation in (|2.?1 ) is elementary in the sense that 



where, as usual, the indicator function 1^ denotes a random variable with value 1 on the event 
A and value on the complement of A. 



3 An algorithm for the VaR contributions 



The aim with this section is to derive expressions for the VaR contributions defined by (2,7) 
which can be evaluated by the usual output of CreditRisk + or at least by minor modifications 
of this output. 



In order to derive an expression for the expectation in (|2.7|) which can be numerically evaluated 



we have to specify the a-algebra from Assumption 2.1 and the stochastic intensities R\, . . . , R n 
from ( |2.2| ). This specification constitutes what is commonly called the CreditRisk + model. Recall 
that the Gamma-distribution with parameters (a, (5) G (0, oo) 2 is defined by the density f a> p 
with 



_a— 1 

x -x//3 



Assumption 3.1 The stochastic intensities R{ of the A, in (B.i ), i = 1, . . . , n, are given as 

k 



R; 



r i' i ^ ' 

3=0 



(3.2) 



where rj t i > 0, J2j=i r j,i > 0' Ya=i r o,i = So > is a constant, and the Si,... ,Sk are 
independent and Gamma- distributed with parameters (ctj,j3j) £ (0, oo) 2 ; j = 1,... , k. The a- 



algebra A is the a -algebra generated by (Si,... ,Sk), i-e. A = a(Si, . . . ,Sj-)- Eq. (2.Z.) holds. 



Usually, the variables Sj are interpreted as systematic default intensities which are characteristic 
for sectors, countries, or branches. Equation ( |3.2; ) then mirrors the fact that a firm may be influ- 
enced by the evolutions in several sectors. Since E[5j] = ay f3j, j = 1, . . . , k, under Assumption 



3A eq. (2.2) is equivalent to 



Pi = r 0>i S + } y r^j ctj i = l,...,n. 



(3.3) 



Remark 3.2 
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(i) In practice, it is unlikely to have a credit portfolio parameterized like in Assumption 3.1 . 



In \CSFB (199% sec. A12) and \Lehrbass et al. (200$ ), for instance, the given data are the 



default probabilities pi, i = 1, . . . , n, and default volatilities o~i, i = 1, . . . ,n, for each credit 
as well as factor weights 9j t i, j = 0, . . . , k, i = 1, . . . ,n which measure to which degree 
credit i is exposed to factor j and satisfy X^j=o = 1> * = 1, • • • > w. /n this case, the 
are defined by 



3,1 



(3.4a) 



with Hj = Yl7=iPi = • Hence, So is given by So = hq. It is quite common to specify 

also the variances tJ of the factors Sj, j = 1, . . . , k in order to determine the values of 
(aj,/3j). Since E[Sj] = ctj f3j and vax[Sj] = otj0j, we have, once the expectations \ij and 
the variances r? are known, 



a ; 



Pi 



j_ 

2 ' 
3 
.2 



(3.4b) 



In CSFB (1991 ), the choice Tj = X^ILi ai ^ s suggested in order to get a link between 
the default volatilities o~i and the factor variances Tj. However, this approach tends to 
underestimate the factor variances (see Kluge and Lehrbass, 2001). Therefore Kluge and\ 



Lehrbast propose to calculate the factor variances by 



'3,1 



Pi, j = 1, 



(3.4c) 



(ii) In contrast to Biirgisser et al. (1999\) and \Gordy ( 200l ), we do not assume ~E[S. 



1 



in (S.i ). Indeed, this would imply ctjj3j = 1 and, as a consequence, would render the 



formulation of the main result in Corollary \3.4 quite difficult. In a self-explanatory mixture 
of our notions and the notions in \Gordy (2001 ), Gordo 's equation (1) and $3.2\ ) above are 
related by 



WijPi 



Ei=i w ijPi 

n 

E i s j] = ^2 w ijPh j = 0,...,k 



, i = 1,... ,n,j = 0, ... ,k, 



(3.5) 



1=1 



From Assumption 34, we obtain the following representation of the generating function g(z) 
E[z L ] of the distribution of L (cf. |Gordy, 2001] [Lehrbass et al., 200 1Q 



g(z) = exp ( S (£r , 4 z n - l) ) J] ( 1 + ft - <£ rj4 

3=1 



i=l 



i=l 



(3.6a) 
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Appealing to the binomial series it is clear that the radius of convergence p g of the power series 
representation of g is given by 



Pg 



sup 



n 1 

{z>0: Y, r ^ zVi < 1 + 73-, J = !,-•• ,^}- 

Pi 



i=l 



(3.6b) 



Obviously, we have p g > 1, i.e. there are numbers z > 1 for which the power series of g converges. 
From (|3.6a ), the exact distribution of the portfolio loss L can be successively determined (see 
CSFB, 1997| , A10) by means of a recurrence relation. 



More important here, we can derive from ( |3.6a ) the generating functions of the sequences I 
E[JVil {L=J} ],i = l,... ,n. 



Theorem 3.3 Let N\, . . . ,N n satisfy Assumption 3.1. Define L by (\S-8j ), g by 6a\ ), and p g 



by ( 3. 61 ). Then the functions /, : (0, p g ) — > R ; i = 1, . . . , n with 



E[Ni 



(3.7a) 



converge and can be represented as 



fi(z) = z Vi g{z)[S [ 



0,1 fei + /3i-AELi^/ 



(3.7b) 



Proof. Note that E[L z l ~\ = Yli^o = I] z l = z ■Jj-(z) is finite for z S (0, p g ) as a consequence 
of general properties of power series. Since we assume i/j > 0, by 



/<(*) = E[N t z L ] < f, rl E[L , 
also the /i(^) are finite for z G (0, p g ). Fix i E {1, . . . , n} and 2; G (0, p g ) and observe that 

d 



dh 



L+hNi\ 



\h=Q 



Nj z n log z. 



(3.8) 



(3.9) 



Recall from section |2| that we assume the v\,... , v n to be positive integers. Nevertheless, ( |3.6a| ) 
holds not only for positive integers u\, . . . ,u n but more generally for any v\, . . . , u n > 0. Thus, 
for small h with \h\ < 



E[z 



L+hNil 



exp(S ( £ r 0) i^+r 0)i ^ +ft -l 



3=1 



i=i,i& 



(3.10) 



By ( |3.10 ), ( |3.9| ) will imply ( 3.7b| ) as soon as it is clear that the order of differentiation and 
expectation in ^■E[z i+ ' liVi ] can be exchanged. But this follows from standard results on differ- 
entiation under the integral (e.g. Durrett, 1996, Theorem A. (9.1)). □ 
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Note that Martin et aL (2001 , translated to our notation) suggest the following approximation 
for the risk contributions to VaR in the sense of ( |2.7| ): 



Pi(L) ps v, 



E[N, exp(s 5 L)] 



(3.11a) 



E[exp(s 5 L)] 

where s$ > 0, the so-called saddle point, is given as the unique solution of the following equation 

E[L exp(s^ L)] 



qs(L), 



(3.11b) 



Efexp^L)] 

with qs(L) standing either for qs(L) or any reasonable estimator of it. Note that the right-hand 
side of (plal) can be expressed as j- ^ log E[exp(z (L + h N{))]\ 

(h=0, z=s$) 

This approach is based on the so-called saddle-point method by means of which quantiles to 
tail probabilities of the distribution of L can be approximated (see Martin et al., 2001 ; GordyJ 



20011) . 



Under Assumption 3.1 , we obtain from Theorem |3.3| a rather simple formula for the right-hand 
side of ( pla| ): 

Pi(L) ps v, expfo s s ) (So r ,i + £ — -= , '' JJ -, A . (3.12) 

V ~[ 1 + Pi ~ Pi z2h=i r j,h exp(^ s s ) ) 



In order to derive the announced algorithm for the calculation of the VaR contributions in (2.7), 
we need a further notation. Note that the probability measure P under consideration depends, 



in particular, on the parameters a\, ... from Assumption 3.1. We express this dependence 
by writing 

P = P a with a = (qi, . . . , ctfc), 
and, analogously, E Q for the expectations. 



Corollary 3.4 Let N\, . . . ,N n satisfy Assumption \3.J\ . Fix a confidence level 5 € (0,1) and 
calculate the value-at-risk qs(L) of the portfolio loss L according to t \2.4\ ) for the probability P a . 
Then, for i = 1, ... ,n, the VaR contributions in the sense of 7p can be calculated by means 
of 



E a [Ni\L = q s (L)] 



S r ,j P a [L = qs(L) - Vi] + Y!j=i a j Pj r j4 F a(j)[ L = Qs{L) - Vi] 



P a [L = q s (L)] 

with a = (ot\, . . . , af.) and a(j) = (a%, ... ,aj + 1, . . . , ctk)- 



(3.13) 



Proof. Fix any i E {1, . . . , re}. By ( |2.8|) , it suffices to show that the numerator of the right-hand 
side of ( p.13 ) equals E^fiVj l{L=q s (L)}]- We will prove the even stronger identity 



Ea[Nil{ L= t}] = S r 0ti P a [L = t - +^aj Pj r j:i P a{j) [L = t - vi] , (3.14) 

3=1 
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for any non-negative integer t. By Theorem we know the generating function fi{z) of the 
sequence 1 1— > ~E a [Ni liL=t}]- Write g a instead of just g in order to express the fact that g like P 
depends on the parameters a\, . . . , a*, from Assumption 3A. With this notation, ( 3.7b|) can be 
written as 



fi(z) = z u * I S r ,i g a (z) + ^ aj 0j r jti g a (j){z) j 
V j=i / 

oo / k \ 

= ^2{ s o ro,i P a [L = t- Uj] + ay /3j r^j P Q(j) [£ = t - J 
t=o ^ 7=1 ' 



(3.15) 



Note that for the second identity in ( p,15| ) we have used the fact that L is non-negative. By the 
uniqueness of power series, ( |3.7a ) and ( |3.15 ) imply now ( 3.14| ). □ 



Remark 3.5 



(i) Corollary 3.4 states essentially that the VaR contributions according to (2.7) can be de- 
termined by calculating the loss distribution (k + 1) times with different parameters. This 



is quite easy if the portfolio is parameterized as in Assumption 3.1. But also in the case 
described in Remark 3.i (i) (with ^3.4-c ) as definition for Tj) it is easy to run the algo- 
rithm in the appropriate parameterizations. Just note that in order to calculate the P a (j) 
probabilities, it suffices to replace the weights 6j t i, i = 1, . . . ,n, with 9'^ 



(ii) By construction of Ni as a conditionally Poisson distributed random variable, we have 
P Q [iVi > 1] > 0. Hence it is possible that the VaR contributions according to ( 3.13Q and 



( 2. % ) become greater than the exposures Vi. 



(Hi) A further consequence of the model construction according to Assumption 3A is that 

E a [Ni\L = q s (L)] = (3.16) 
may happen. In fact, for fixed i E {1, ... ,n} the following statement holds: 

n 

E Q [iVj|L = qs(L)] equals if and only if Yl u j u j = Qs{L) with non-negative integers 

i i 

u\, . . . ,u n implies Uj = 0. 

To see this, define Y = L — z/j Ni, and observe that 



P[L = q 5 (L) -Vi]=0 



P[Y = q s (L) —kui] = 0, 1 < k < q s (L)/ui 
P[L = q 5 (L)] = P[Y = q 6 (L),Ni = 0] 
^ {L = q 5 (L)} = {Y = q s (L), N t = 0}. (3.17) 

Denote by M the set of all n-tuples (ui, . . . ,u n ) of non-negative integers such that 
5^j=i u j u j = Qs(L). Then the assertion follows from ( \3.1 r / ) and the fact that 

{L = q s (L)} = II {jVi =«i,... ,N n = u n }. (3.18) 



U {^1 

(ttl,... ,U n )GM 



In particular, ft3.lt ) can occur in case > qs(L). This is possible whenever pi < 1 — 5 
and vi is large compared to the other exposures. 
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Remark ^ (ii) is being caused by the fact that CreditRisk" 1 " makes use of a Poisson approxi- 
mation in order to compute the portfolio VaR. Hence one has to consider the Poisson variables 
Ni when computing risk contributions. As a consequence, there is little hope that switching to 
another risk measure will solve the problem of contributions which are larger than the corre- 
sponding exposures. The situation concerning Remark |3.5| (iii) is more favorable. 



For instance, one can try to smooth the conditional expectations in ( 3.13 ) by choosing a fixed 
positive integer t and computing E Q [jVj | qs(L) — t<L< q$(L) + t] instead of E a [Ni \ L = qg(L)]. 
We will not go into the details of this approach since in the subsequent section we will see, 
that, in particular, the problem from Remark |3.5| (iii) can be avoided by switching from VaR 
to Expected Shortfall, a risk measure which should be preferred for some reasons ( |Acerbi and 



Tasche, 2001). 



4 Contributions to Expected Shortfall 

The Expected Shortfall (ES) at level 5 G (0, 1) of the portfolio loss L can be defined as 

ES 5 (L) = (1-8)- 1 j q u (L)du. (4.1) 

It may be characterized as the smallest coherent risk measure dominating VaR and only de- 
pending on L through its distribution ( pelbaen, 199S , Th. 6.10). Of course, (fO|) appears not 



very handy for calculations or defining ES contributions. Nevertheless, if L had a continuous 
distribution, the representation 

ES S (L) = E[L\L>q 5 {L)} (4.2) 



would be equivalent to (4.1) ( Acerbi and Tasche, 200i~| , Cor. 5.3). Indeed, if we assume to deal 



with a large portfolio we can hope that the difference between ( [4.1[) and (|4.2| ) will not be too 



large. With (4^2) as definition, the decomposition 

n 

E[L\L>q s (L)} = Y,Vi E [ N i\ L ><ls(L)) (4-3) 



suggests the choice of Vi E[iVj | L > q$(L)] as ES contribution of asset i. Further evidence for this 
choice stems from the fact under certain smoothness assumptions ( |4.3j ) is just the decomposition 
corresponding to (^) with p = ES flTasche, 199£ , Lemma 5.6). 



From Corollary we obtain the following result on the computation of the ES contributions. 
Corollary 4.1 Let N%, . . . ,N n satisfy Assumption [O. Fix a confidence level 5 G (0,1) and 



calculate the value-at-risk qs(L) of the portfolio loss L according to (2.4) for the probability P a . 



Then, for i = 1, . . . ,ra, the ES contributions in the sense of can be calculated by means of 

„ , AT , T . S o r o,i P*[L > qs(L) - Vi] + £ fe =1 Qj (3 j r jti P a(j) [L > q s (L) - Vi ] 

E a[Nl | L > q s (L)] = P a [L>q s (L)} ' ^ 

with a = (ai, . . . , «fc) and a(j) = (a%, ... , ctj + 1, . . . , at)- 
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E^M | L > , 5 (L)] = Pa[L > g , (L)] ( 4 - 5 ) 



Proof. Observe that 



and note that in the proof of Corollary p.4| we have in fact shown that 

k 

E a [Ni l{L=t}] = S r , i P a [L = t-u i } + ^2a j P j r jii P a{j) [L = t-u i } (4.6) 

i=i 

for any non-negative integer t. □ 



Note that none of the probabilities in the numerator of the right-hand side of ( |4.4j ) can become 
0, since by definition of q$(L) the probability P a [L > q$(L) — Vj\ satisfies 

V a [L>q s {L)- Vi \ > P a [L>q s (L)] > 1-6 (4.7a) 

and for any integer t we have 

P a [L = t]=0 P a{j) [L = t) = 0, j = 1, . . . , k. (4.7b) 

Nonetheless, the observation in Remark |3.5| (ii) holds for the ES contributions as well. 



5 Conclusion 

Attributing the total risk of a credit portfolio to its components in an exhaustive and risk 
respecting way is an important task in portfolio management. For the case that risk is measured 
as value-at-risk (VaR) and determined with the CreditRisk + model ( CSFB, 1997| ), we have 



shown that this attribution can be performed by calculating the loss distribution {k + 1) times 
{k denoting the number of sectors in the model) with slightly different parameters. The same 
statement holds for the attribution when risk is measured as Expected Shortfall (ES). 
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